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The relation between rotating plane Couette and Taylor-Couette flows is clarified. Experimental 
data are used to quantify the behavior of the minimum Reynolds number for subcritical turbulence 
as a function of rotation and curvature. This last dependence is understood through a phenomeno- 
logical analysis, which also allows us to relate the subcritical turbulent transport efficiency to the 
transition Reynolds number. This implies that the Coriolis force reduces the efficiency of subcritical 
turbulent transport with respect to nonrotating flows, and resolves an ongoing controversy. 
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Shear flows constitute one of the prototypical type of 
hydrodynamical flows. Furthermore, they are commonly 
found in various instances, e.g. geophysical and astro- 
physical contexts. This makes the understanding of their 
properties an important issue from both a fundamental 
and a practical point of view. Most prominently, char- 
acterizing turbulence and turbulent transport in such 
flows is critically needed, as astrophysical and geophysi- 
cal shear flows are usually fully turbulent because of the 
large scales involved. 

A large body of experimental evidence has been col- 
lected on Taylor-Couette flows in the linearly unstable 
regime (see, e.g., Ref. Q and references therein). Much 
less is known about the linearly stable, or subcritical, 
regime 0, 01 Plane Couette flows (rotating or not) 
are more difficult to realize experimentally, and have been 
less extensively studied 0, 0, • 

On the theoretical side, intense efforts have been de- 
voted to the understanding of turbulence spectral and 
statistical properties (see, e.g., || for an introduction to 
the subject); also, mostly for practical purposes, rather 
complex turbulent transport models have been developed 
0. Unfortunately, in spite of these remarkable successes, 
some basic properties such as magnitude of the critical 
Reynolds number of fully developed turbulence, or its 
dependence on rotation, are not yet understood. Worse, 
even the most sophisticated Reynolds stress closure mod- 
els fail to account for the existence of subcritical turbu- 
lence in the presence of a stabilizing rotation. 

Some of these shortcomings are addressed in the 
present investigation, mostly through a phcnomcnolog- 
ical analysis of the effects of rotation and curvature in 
rotating plane Couette flows and Taylor-Couette flows. 
This provides us with an understanding of some of the 
characteristic features of subcritical shear turbulence, 
and, most importantly from a practical point of view, es- 
tablishes a relation between the efficiency of subcritical 
turbulent transport and the magnitude of the minimum 
turbulent Reynolds number. The emphasis on subcriti- 
cal turbulence follows for the following reasons. First, an 
analysis of angular momentum transport suggests that 
both linearly stable and unstable fully turbulent Taylor- 



Couette flows are controlled by similar nonlinear physics 
0, E| ■ Secondly, it turns out that subcritical shear flows 
are easier to analyze from a phenomenological point of 
view, and they are in any case important in themselves. 

Let us first reexamine the connection between rotat- 
ing plane Couette flows and Taylor-Couette flows. The 
Navier-Stokes equation for rotating plane Couette flows 
reads, in the rotating frame 



dw _ 

— h w. vw = 

ot 



VP* 



— 2J~i x w + vAw. 



(1) 



In this equation, P* represents the sum of the inertial 
term and of the gas pressure term, as usual when con- 
sidering incompressible flows. Such flows are character- 
ized by two dimensionless numbers, the Reynolds number 
Re = AVAL/v ~ \w.Vw /vAw\, and a rotation (inverse 
Rossby) number R n = e2flAL/AV ~ \2fl X w/w.Vw\, 
which is a global measure of the rotation parameter de- 
fined by S = —2fl/(dV/dy). In these definitions, the x 
axis is assumed to lie in the streamwise direction, the 
y axis in the shearwise direction, and the z axis in the 
spanwise direction (rotation axis); V(y) is the mean flow 
(along x). The rotation number is positive (resp. neg- 
ative) (e = ±1) depending on whether the rotation is 
cyclonic (resp. anticyclonic) . 

Stability limits for inviscid rotating plane Couette 
flows can be determined through a displaced particle 
analysis 0, 0] . Indeed, although the total work of the 
Coriolis force vanishes, the work of the Coriolis force com- 
ponent in the streamwise direction during a fluid particle 
displacement in the shearwise direction does not. Com- 
paring the resulting post-displacement streamwise veloc- 
ity to the equilibrium one implies that the flow is linearly 
unstable if — 1 < S < 0, and linearly stable otherwise, as 
the squared oscillation frequency of fluid particles is given 
by (dV/dy) 2 S(S + l). 

The Navier-Stokes equation for Taylor-Couette flows 
is most meaningfully compared to that of rotating plane 
Couette flows when it is dispelled in a frame rotating 
with some characteristic rotation velocity of the flow fl a 
(e.g., the average angular velocity of the two cylinders), 
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as only differential rotation plays a role in the generation 
of turbulence. Designating by w and <j> the velocity and 
azimuthal coordinate in the rotating frame, the Navier- 
Stokes equation for w — (w r , «u, w z ) becomes 



duo 

-— + w.V'w + 2f2 Q x w- 

ot 

w% 2 W(t> w r VP* 

— -e r H - — ej, = h uAw, 

r r p 



(2) 



where P* includes again all gradient terms, and 



w.'V'w = (w.Vio r )c r + [rw.'V{w ( f 1 /r)\e ( f ) + 

(w.Vw z )e z . (3) 

With this new definition, the contribution of the az- 
imuthal velocity derivative to the "advection" term van- 
ishes when the fluid is not sheared; as a consequence, 
the w. V'w can appropriately be named the "advection- 
shear" term. This reexpression of the Navier-Stokes 
equation is essential: explicitly isolating shear terms in 
the advection term distinguishes the effect of the shear, 
upon which the Reynolds number is built, from the oth- 
ers, the flow curvature and rotation. Indeed, with respect 
to rotating plane Couette flows, Eq. QJ, two new terms 
are added, Wg/r and 2w r wg/r, which represent the ef- 
fect of the curvature of the flow due to the cylindrical 
geometry. As will soon become apparent, Eq. © is ex- 
tremely helpful to develop a correct physical analysis of 
Taylor- Couette flows. 

This reexpression of the Navier-Stokes equation shows 
that the three dimensionless numbers which best char- 
acterize Taylor- Couette flows are the Reynolds num- 
ber Re = RAVtAR/v - \w.V'w/vAw\, the rotation 
number R n = e2VtAR/RAVt ~ \2Q, x w/w.V'w\, 
which is a global measure of the rotation parameter 
S = 2Q I (rdil / dr) |23j, and the curvature number R c = 
AR/R ~ w 2 /(r\w. V'u>|, which is a global measure of 
the curvature parameter C — 2wg/(r 2 d£l/dr). As will be- 
come apparent shortly, these numbers are to be preferred 
over the most widely used inner and outer Reynolds num- 
ber, and ratio of cylinder radii rj. 

Note also that rotating plane Couette flows can be 
viewed as a limiting case of Taylor-Couette flows, when 
R c — > at finite Rq. This implies that there must 
be a connection between the stability limit of rotating 
plane Couette flows recalled above, and the Rayleigh 
criterion. Indeed, the displaced particle analysis pre- 
sented above is easily transposed to Taylor-Couette flows 
[l3|. Let us consider a fluid particle at some arbitrary 
radius r a , and chose Q = 0,(r o ). With this choice, 
the curvature terms are second order with respect to 
the radial displacement of the particle, and can be ne- 
glected. Therefore, the argument for rotating plane Cou- 



ette flows applies in exactly the same way to Taylor- 
Couette flows, and implies that they are linearly unstable 
when — 1 < S = 20 / (rdCl / dr) < 0, where S is evaluated 
at r = r a . Equivalently, the restoring (or destabilizing) 
frequency reduces to the well-known epicyclic frequency: 
{rdn/dr) 2 S(S+l) = {2fl/r)d(r 2 n)/dr. This implies that 
stability follows when the Rayleigh criterion is satisfied, 
as expected, and shows that this criterion can be reinter- 
preted in terms of the action of the Coriolis force in the 
rotating frame. 

A similar relation exists between the two types of flows 
in what concerns the minimum Reynolds number for the 
onset of subcritical turbulence. However, this quantity 
for Taylor-Couette flows depends both on rotation num- 
ber and curvature number, an important feature whose 
meaning and consequences have apparently escaped at- 
tention until now. 

Transition to turbulence in subcritical flows occurs 
abruptly, without undergoing the series of bifurcations 
which is characteristic of linearly unstable flows. Note 
however that the minimum turbulent Reynolds number 
is not defined in a clear-cut way, as the transition from 
laminar to turbulent occurs in general at higher Reynolds 
number than the transition from turbulent to laminar, 
and as there is a range of Reynolds numbers over which 
the flow changes from highly intermittent to fully turbu- 
lent. Nevertheless, the data quoted here are more or less 
directly comparable, as they are rather characteristic of 
fully turbulent flows, except for the data of Ref. [4[ , which 
are therefore used in a qualitative rather than quantita- 
tive way in the discussion, when needed. 

The available data on cyclonically rotating plane Cou- 
ette flows imply that the minimum turbulent Reynolds 
number increases with the rotation number as 61 



Re m ~ Re c + a c R n , 



(4) 



where Re c ~ 1400 is the transition Reynolds number of 
plane Couette flows ja], and a c ~ 26000. The data were 
obtained for low cyclonic (0.1 > Rq) rotation. 

The experimental results of Wendt 0] and Taylor 
for Taylor-Couette flows in cyclonic rotation were ob- 
tained by maintaining the inner cylinder at rest. In this 
case, R c = Rn — AR/R, which makes it a priori diffi- 
cult to distinguish the effects of rotation from those of 
curvature. Nevertheless, for AR/R < 1/20, the data 
are clearly consistent with Eq. (0J (in agreement with 
the existence of the rotating plane Couette flow limit for 
Taylor-Couette flows), while for AR/R > 1/20, the min- 
imum Reynolds number varies as Re m ~ Re* (AR/R) 2 , 
with Re* ~ 6 x 10 5 01- This last behavior is necessarily 
due to the effect of the flow curvature, and not rota- 
tion. This follows first because the predictions of Eq. Q 
fall way below the required minimum turbulent Reynolds 
number when extrapolated to larger AR/R, and because 
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the experiments of Wendt for marginally stable anticy- 
clonic flows are globally consistent with the quadratic 
dependence of Re„. for large enough AR/R 13]; also, 
the results of Ref. |j| are in qualitative agreement with 
this conclusion All these arguments imply that, for 
cyclonically rotating Taylor-Couette flows, 



Re r 



Re 



c Rn + Re*R 2 . 



(5) 



This relation applies to moderate rotation rates (Ro < 
1/2), where the extrapolation of Eq. (J3J is still valid 
and the turbulence still three-dimensional. The fact that 
the rotation and curvature contributions add is implied 
by the data of Richard Q (see Ref. for details). 
Notice that the contributions of either rotation (a c Rn) 
or curvature {R*R 2 ) become comparable to Re c when 
the rotation or curvature number is ~ 1/20. This re- 
markable feature is no coincidence. Indeed, the tur- 
bulence regeneration mechanism identified in nonrotat- 
ing plane Couette flows 0, 0] has an overall timescale 
t turb ~ \{)Q{dU / dy)' 1 ~ lOOirdn/dr)- 1 . The timescale 
associated with the Coriolis term is in ~ f2, while the 
timescale associated with the curvature terms is t c urv ~ 
r/wg ~ Ai7 _1 . Therefore, these effects are expected to 
significantly affect the turbulence regeneration mecha- 
nism when these timescales decrease to become compa- 
rable to a few times t tur b, or equivalently, when Ro or R c 
exceed a few percents. This physical constraint is what 
primarily determines the magnitude of a c and Re* (once 
the form of the dependence on i?n and R c is known) . In- 
deed, for R c , Rq > 1/20, requiring a c Rn, Re*R 2 > Re c 
implies that a c ~ 10 4 and Re* ~ 10 5 . 

The explanation of the dependence of Re m on R r of 
Zeldovich 0] is inconsistent with some of the data [l(| , 
while the explanation of Dubrulle 01 is incompatible 
with the fact that this is an effect of the flow curva- 
ture. In fact, it can be understood by extending the 
turbulent viscosity description 0], an argument which 
is clarified and noticeably improved here. Characterizing 
turbulent fluctuations responsible for turbulent transport 
by their coherence scale Im and velocity fluctuation vm, 
a simple analogy with molecular transport implies that 
[lf| {SvrSv^) ~ v t rd(Vl) /dr, with v t ~ Imvm- In a Kol- 
mogorov cascade picture, the rate of energy dissipation 
e can be expressed in terms of the same quantities as 
e ~ C*v 3 M /l M ^ (is t /2)(rd(n}/dr) 2 , where C* ~ 0.1 
is a rather universal constant In this picture, the 
turbulence regeneration mechanism involves scales in the 
range AL to Im, which do not dominate the turbulent 
transport; furthermore, anisotropy is neglected for the 
moderate shear and rotation numbers involved. Finally, 
a flow submitted to a velocity shear tries to restore global 
thermodynamic equilibrium by suppressing the shear and 
transporting momentum in the shcarwise direction. A 
linearly stable system has only two means to achieve this 
purpose, laminar or subcritical turbulent transport, and 



will chose the most efficient one. Turbulent transport is 
dominant when \{5v r 5v,f,)\ > i/\rdfl/dr\, so that, at the 
transition to turbulence, v t — C v v where C„ is a constant 
of order unity. Taken together, these relations imply that 



VM 



(2C*) 1 /2 



1/2 



rAfl 



1/2 ' 



l M * (c^C*) 1 ' 2 



Re 



1/2 Ar 



Re 



1/2 ' 



(G) 



(7) 



where Ar is the typical scale of the shear away from the 
boundaries, and rAQ the corresponding shear amplitude. 

Let us focus for the time being on Taylor-Couette sys- 
tems dominated by the curvature. By making the outer 
cylinder radius arbitrarily large, this last relation implies 
that, if Re m were constant (independent of the relative 
gap width), the turbulent scales could become arbitrar- 
ily larger than the radius, at any given radial location 
in the flow, which makes no sense. Reversely, this shows 
that the local radius r, at any flow location, must impose 
a limiting scale. Indeed, by imposing Im — 7?" (where 
the proportionality constant 7 will soon be determined) , 
Eq. is equivalent to [24| 



Re, 



C^C*) 1 ' 2 (ArY „ . ( Ar 



r 



- = Re " I — I ■ (8) 



In the same vein, the resulting turbulent velocity ex- 
pression vt ~ {C v j ' Re*)r 3 \Att/ Ar| = (3r 3 dil/dr also im- 
plies that in this regime, turbulence is controlled by local 
rather than global conditions in the bulk of the flow. This 
turbulent viscosity prescription has indeed been derived 
from the data of Wendt [l(|, which imply that ~ 10~ 5 , 
so that C v ~ 6 and 7 ~ 1/500. 

If turbulence in the bulk of the flow is controlled by 
local conditions, this implies that the width of the flow 
plays no role in the curvature dominated regime (except 
indirectly through the way global conditions might con- 
trol local ones). This surprising conclusion can be under- 
stood in the following way. For narrow gap widths (still 
neglecting the effect of rotation), the flow reduces to a 
(nonrotating) plane Couette flow, and turbulence occurs 
as soon as RAQAR/v exceeds Re c . When AR > R/20, 
the mean radius R locally defines radial "boxes" of size 
i?/20 which are turbulent if (r dtt / dr) {R / '20) 2 jv > Re c , 
so that, for the shear to be large enough on scale R/20, 
the Reynolds number must increase quadratically with 
the gap width. Incidentally, this shows that the tur- 
bulence regeneration mechanism in curvature dominated 
flows must be closely related to the one that operates 
in nonrotating plane Couette flows. Consequently, the 
mechanism identified in Ref. [Til Il5| is probably more 
than a near- wall turbulence regeneration mechanism. 
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Note in this respect that this mechanism is largely in- 
sensitive to the actual nature of the boundary condition; 
only a box size is needed. 

What has been achieved so far ? The magnitude of Re c 
is explained by the regeneration mechanism just men- 
tioned, while the magnitudes of a c and Re* have been 
constrained from the overall time scale of this mecha- 
nism; the meaning of the quadratic dependence on R c 
has also been brought to light. Only the linear depen- 
dence of Re m on Rq still needs to be explained. As the 
Coriolis force term does not pinpoint a particular length 
scale (only a time scale is involved), in contrast to the 
curvature term, the flow rotation is expected to affect 
the regeneration mechanism in an important way, so that 
this behavior is beyond the scope of the phenomenolog- 
ical analysis developed here. It is essential to note that 
the turbulent transport v t ~ ImVm based on Eqs. © and 
Q is valid for all Re ^ Re m , as long as the turbulence 
regeneration mechanism does not change with increas- 
ing Reynolds number; this follows because increasing the 
Reynolds number only increases the range of the inertial 
spectrum, while the regeneration mechanism necessarily 
acts on larger scales 0, This feature always al- 

lows us to relate the minimum Reynolds number to the 
transport efficiency through Eqs. J^J and Q and the tur- 
bulent viscosity prescription, an essential feature from a 
practical point of view. 

Finally, the dependence of Re m on Rq allows us to 
resolve an important controversy on the relevance of 
subcritical shear turbulence to accretion disk transport. 
Most astrophysicists now believe that the Coriolis force 
relaminarizes shear flows in the anticyclonic regime (in 
particular in the "keplerian", Rq = —4/3, regime), on 
the basis of two sets of simulations [20L l2l| performed 
for rotating plane Couette flows with pseudo-periodic 
boundary conditions and with a tidal term included in 
the generalized pressure P* . On the other hand, a mi- 
nority insists that this point of view is inconsistent with 
the experimental evidence 01 ■ This last claim is indeed 
correct 0,E1; nevertheless, this can be reconciled with 
the simulation results as is now argued with a much more 
focused argument than the one given in Ref. |18| . First, 
although no usable data are available in the relevant ro- 
tation dominated anticyclonic regime, the approximate 
"Richardson similarity" of rotating flows |22j (valid for 
the low enough rotation numbers of interest here) im- 
plies that the minimum Reynolds number of anticyclonic 
flows (Rq ^ — 1) behaves as Re m ~ Re c + a ac |l + Rq\ 
with a ac > a c (see Ref. [l3j for details). This im- 
plies that the turbulent amplitude decreases steeply with 
the rotation number, as implied b y th e simulations per- 
formed close to marginal stability |2(| . Furthermore, for 
Rq = -4/3, AR/l M ~ (AV)/v M > a few 100, and 
(SvrSvrf,) ~ 10~ 3 (flAR) 2 . These figures imply that the 
largest simulations of keplerian flows performed to date 



21] miss the mark by lack of resolution, but not by much. 

The numerical results were most probably not recog- 
nized as characteristic of turbulence, because it is widely 
believed that turbulent scales and velocity fluctuations 
are comparable to the scale of the flow, and mean velocity 
difference on this scale. On the contrary, the analysis pre- 
sented here implies that rotation reduces the transport 
generating velocity fluctuations with respect to nonrotat- 
ing flows. A fuller account of the relevance of these fea- 
tures to accretion disk turbulent transport will be given 
elsewhere 1 13. 
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